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ABSTRACT 

Most studies of Non-Gaussianity (NG) in the CMB data rely on moment based approaches that depend 
on the study of bispectrum or its higher order analogs of multispectra. In contrast, the studies that use 
the Minkowski functionals (MF) depend on morphological characteristics of the observed fields. The 
two approaches are complementary. Indeed, in a perturbative expansion of MFs, the lowest order terms 
correspond to the angular bispectrum albeit with a different set of weights for the individual modes. For 
all angular scales study of MFs are equivalent study of three one-point generalized skewness parameters. 
There has been lot of recent progress in estimation of skewness mainly through estimation of a spectrum 
associated with the ordinary skewness (namely the skew-spectrum). We show how the topological statis- 
tics such as MFs can also be computed from real data sets in the presence of mask and inhomogeneous 
noise by using methods for estimation of skew-spectrum. We exploit the fact that the computation of 
MF at the lowest order is equivalent to the estimation of three different skewness parameters. These one- 
point estimates are volume averages of third order statistics of fields that are constructed from the original 
data. Generalizing the concept of ordinary skew-spectra, which was recently introduced for the study of 
bispectrum, we define a sets of three different skew-spectra which specify the MFs at lowest order in 
non-Gaussianity and can probe them as a function of harmonic number. These spectra can also be studied 
independently to provide a valuable check for any cross-contamination from the secondaries or fore- 
ground sources. The Pseudo-C; (PCL) approach is employed to estimate the generalized skew-spectra 
associated with the MFs in the presence of mask from noisy data. The variance for such estimators is 
analyzed. The results presented here are generic and can be useful in analyzing data from other projected 
surveys such as from the weak lensing surveys or from the future Sunyaev-Zel'dovich (SZ) surveys. Gen- 
eralization to 3D can be done in a straight forward way which will be useful in quantifying the topology 
of galaxy distributions. We also go beyond the lowest order in skew-spectra and discuss the extraction of 
four generalized kurtosis and the corresponding kurt-spectra that are relevant in studying the next order 
corrections. 
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1 INTRODUCTION 

The study of Cosmic Microwave Background (CMB) radiation provides the cleanest window to probe the very early stages of the Universe's history. This 
can be used to probe the mechanism that generates seed perturbations which leads to the structure that we observe in the present day Universe. Recent 
observations by WMAP0 supports the basic predictions of inflationary scenarios. It confirms that the primordial perturbations are nearly scale-invaria nt, 
adiabatic and nearly Gaussian. It is expected that the ongoing PlanclQ surveyor will improve our knowledge in these areas (see e.g lTauber etHI bold) ). 

It is well established now that non-Gaussianity from simplest inflationary models based on a sing le slowly-rolling s calar field is typically very smal l 
dSalopek & Bondlll99(]|.[l99lMFalk et al.ll993klGangui et al.ll994l ; lAcauaviva et al.l200l lMaldacena 20031) . (see e.g. lBartolo. Matarrese & Riottol (120060 



for a review). However there are man y variants of simple inflationary models which include models with multiple scalar fields ( Linde & Mukhanov 



1 19971; iLvth. Ungarelli & Wandsl)2003l). features in the inflationary potential, non-adiabatic fluctuations, non-standard kinetic terms, warm inflation 

I0J; 



( Gupta . Berera & Heavens! I2OO2I : iMoss & Xiond l2007t) . or deviations from Bunch-Davies vacuum that can all lead to much higher level of non- 
Gaussianity. 

Early observational work on the bispectrum from COBE dKomatsu et alj|2002h and MA XIMA dSantos et alj |2003) was followed by much more 



accurate analysis with WMAP dKomatsu et alJl2003l ; ICreminelli et al.ll2007 : Spergel et al. 2007). With the recent claim of a detection of non-Gaussianity 



1 http://map.gsfc.nasa.gov/ 

2 http://www.rssd.esa. int/index.php?project=Planck 
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jYadav & Wandeltll2008h in the Wilkinson Microwave Anisotropy Probe 5-year (WMAP5) sky maps, interest in non-Gaussianity has obtained a tremen- 
dous boost. Much of the interest in primordial non-Gaussianity h as focused on a ph enomenological 'local /jvl' parametrization in terms of the perturba- 
tive non-linear coupling in the primordial curvature perturbation d Verde et all2007t) : 

= * t (a0 + " <*£(*)» + 9NL$l(x) + h NL {$l{x) ~ 3<$|(:r)> 2 ) + ■ • • , (1) 

where <&l[x) denotes the linear Gaussian part of the Bardeen curvature and /jvl, 9nl, Jinl are the non-linear coupling parameters. A number of 
models have non-Gaussianity which can be approximated by this form. The leading order non-Gaussianity therefore is at the level of the bispectrum, 
or in configuration space at the three-point level . Many studie s involving primordial non-Gaussi anity have used the bispectrum, mo t ivated by the fact 
that it contains all the information about Jnl d Babicb 20 051). It has been extensively studied jKomatsu. Spergel & Wandeltll2"005l ; ICreminellil 120031 ; 
ICreminelli et~al]|2006l ; iMedeiros & Contaldill2006l : ICabella et al. 1120061 ; lliguori et all [20071; ISmith. Senatore & Zaldarriagall2009h . with most of these 
measurements providi ng convolved estimates of the bispectrum. Optimized 3-point estimators were introduced b y Heavens ( 1998), a nd have been suc- 
cessively developed jKomatsu. Spergel & Wandeltll2005l ; ICreminelli et alj|2006l : ICreminelli, Senatore. & Zaldarriagal l200l ISmith Zahn & Dorell2000l : 
ISmith & Zaldarriagal 120061) to the point where an estim ator for /jvl which saturates the Cramer-Rao bound exists for partial sky coverage and inho- 
mogeneous noise l lSmith. Senatore & Zaldarriagall2009l). Approximate forms also exis t for equilateral non-Gaussi anity, which may arise in models with 
non-minimal Lagrangian with higher-derivative terms dChen. Huang & Kachrull200^ ; IChen. Easther & Limll2007t) . In these models, the largest signal 
comes from spherical harmonic modes with l\ ~ £2 — £3, whereas for the local model, the signal is highest when one £ is much smaller than the other 
two - the so-called squeezed configuration. 

W hile the bispectrum or its higher order analogs, multispectra, are more commonly used in studying departure from Gaussianity dBartolo et al.l 
2004), alternative statistic s such as Minkowski functional (MFs) too are routinely used for this purpose. MFs describe th e morpholog ical features of a 
fluctuating (random) fields lMecke, Buchert & WagneJ l fl994l) : ISchmalzing & Bucher3(ll997l);ISchmalzing & Gofskil l fl99ll) ; IWinitizki. Kosowskvl ( ll998h . 
The MFs for a Gaussian field is well understood and close form expressions exist jTomitall 19861) . Any departure from the MFs of a Gaussian field wi ll 
signify detection of non-Gaussian ity. Th e MFs have been used to detect non-Gaussianity using projected (2D) fields such as CMB dHikage et alj|2 006). 
Weak Lensi ng (Matsubara & Jai nll200ll) a s well as using 3D d ensity fields as mapped by the Galaxy surveys e.g. SDSS jPark et al.ll2005l ; iHikage et al.l 
I2008L F2006 ; Hik age. Taruva & Sutoll2003l ; IHikage et alj|2002l) . Sev eral analytical resu lts exist for prescriptions to model non-linear gravity as well as 
biasing schemes bo th in quasilinear an d highly nonlinear regimes iHikage et al.l 20081). The MFs have been us ed also for the study of CMB data e.g. 
BOOMERanG data lNatoli eTal] feOld) as well as for the WMAP data jKomatsu et aLKOO^lEriksen et alj|2004l) . The MF based approach has also been 
used to study the effect of lensing on the CMB d Schmalzing. Takada & Futamase 200(Jh 

While the study of primordial NG may be the primary motivation behind t he study of MFs in the context of CMB , they have also been applied to 
probe gravity induced secondary NG using weak-lensing convergence or K-maps Matsubara & Jainl d200ll) ; lTaruva et alj d2002l) . 

One of the main motivations behind studying various alternatives that probe primordial non-Gaussianity has to do with issues related to estimation. 
Different probes are affected differently by different contamination such as presence of foreground or say secondary non-Gaussianity. The methods based 
on multispectra that are typically employed for the estimation of non-Gaussianity typically employ a Fourier (or harmonic) space based approach. On 
the other hand, the techniques developed for estimati on of MFs are tr ad itionally employe d in real space. Previously there has been many attempts to 
probe primordial non-Gaussianity using MFs (see e.g. iNovikovl feOOCh ). jMatsubarall2003h obtained a closed form expression for MFs in D dimension 
using a preservative expansion in terms of various orders of multispectra. At the lowest order, a departure from Gaussianity is characterized by three 
different skewness parameters S (0 \ S (1 \ S 1 ' 2 '. The skewness parameter 5^°' is the ordinary skewness that is most commonly used in various studies of 
non-Gaussianity for projected surveys as well as in 3D. The other set of skewness parameters are defined in terms of different cubic statistics that are 
constructed from the original data using differential operators. 

This paper is organized as follows. In §2 we review the formalism of Minkowski Functionals. In §3 we focus on the one-point estimators, the 
generalized skewness and their links to Minkowski Functionals. In §3 we also introduce the concept of generalized skew-spectra and show how a power 
spectra can be defined for each of the Minkowski Functionals. In §4 the corrections to the lowest order in non-Gaussianity are discussed. These corrections 
are related to estimation of trispectra. In §5 real-life issues such as the mask and noise are discussed and estimators are designed for estimation of the 
Minkowski Functionals using a PCL based approach. In §6 we specialize the general discussion for the CMB sky. Various early Universe models and 
their predictions for the lower order multisp ectra are presente d. Finally, §7 is left for discussion. 

Throughout in this paper a WMAP7 dLarson et alj|20ld) background cosmology will be used for computation of various spectra. The values of the 
parameters /jf£ and /J^ 1 , to be defined later, are set to unity. Some of the results developed here will have wider applicability, e.g. in the context of 
weak lensing surveys or projected galaxy surveys. Detailed results will be presented elsewhere. Results of implementation of our method on WMAP7 
will be presented in an accompanying paper. The techniques developed here can also be applied for the study of SZ maps from future surveys. 



2 FORMALISM 

The MFs are well known morphological des criptors which a re used in the study of random fields. Morphological properties are the properties that remain 
invariant under rotation and translation (see Hadwiger ( 1959) for more formal introduction). They are defined over an excursion set E for a given threshold 
v. The three MFs that are defined for two dimensional (2D) studies can be expressed as: 

V (u) = f da; Vi(y) = i / dl; V 2 (u) = 1. / ndl (2) 
is 4 J as 271 J as 

Here da, dl are the elements for the excursion set E and its boundary <9E. The MFs Vk(v) correspond to the area of the excursion set E, the length of its 
boundary 9E as well as the integral curvature along its boundary which is related to the genus g and hence the Euler characteristics %■ 

In our analysis we will consider a smoothed random field ^(fi) with mean (\I/(f2)} = and variance <7q = (* 2 (f2)}.Fora generic 2D weakly non- 
Gaussian random field ^ on the surface of the sky. The spherical harmonic decomposition using Yi m (Cl) as a basis function = Xw m ^ irnYi m (Q,) 
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Figure 1. In the left and middle panel we display the One point third order statistics, i.e. the three skewness parameters S l (9 B ) as a function of the beam FWHM 8 B . The 
left panel correspond to the local model of primordial non-Gaussianity and the middle panel depicts the equilateral model. In the right panel, the second order one-point 
parameters uq and <j\ as defined in Eq.@ are plotted as a function of the beam FWHM. A WMAP7 background cosmology is assumed. Various plots in each panel 
correspond to the different FWHM as indicated. We consider four different beams with FWHM 9 S = 10', 25', 50', 100'. Note that, being second order in nature, the 
parameters ai do not depend on the modeling of the bispectrum, which is a third order statistics. They depend only on the beam fej and the underlying power spectra C;. 
They are used in normalization of the skew- and kurt-spectra. 




Figure 2. The contributions to the second order one-point parameters [<Tq]j and [uf]i as defined in Eq.@ are plotted as a function of the harmonics /; here = 
X];(2i + ■ A WMAP7 background cosmology is assumed for this calculation. Various plots represent different FWHM as indicated. We consider four different 

beams 8 S = 10', 25', 50', 100'. Being two-point objects, these parameters do not depend on the modeling of the bispectrum and depend only on the experimental beam b; 
and the underlying power spectra C; . They are used in normalization while constructing the three different skew-spectra. 



can be used to define the power spectrum Ci which is sufficient characterization of a Gaussian field im^ i> m ') = Ci^u'^mm' ■ For a non-Gaussian 
field the higher order statistics such as bi- or tri-spectrum can describe the resulting mode-mode coupling. Alternatively the topological measures such 
as Minkowski functionals which include Euler characteristics or genus can be employed to quantify deviation from Gaussianity. Indeed it can be shown 
that the information content in both descriptions are equivalent. At the leading order the MFs can be constructed completely from the knowledge of the 
bispectrum alone. We will be studying the MFs defined over the surface of the celestial sphere but equivalent results can be obtained in 3D using F ourier 
decomposition (Munshi 2010, in preparation). The MFs for a random Gaussian field is well known and is given by Tomita's formula (Tomita 198oj). The 
notations and analytical results in this section are being kept generic however they will be specialized to the case of CMB sky in subsequent discussions. 
The MFs denoted as Vk(y) for a threshold v — \&/<jo; where <tq = (*^ 2 } are defined as follows: 



V k {v) = 



UJ2 



( 27r )(fc + D/2 o; 2 „ fetJfe 



exp - — 



s/2< 



(Jo 



H k {v) 



k(k-l) s(2)Hk 
6 



°> = h ^ {21 + mi + 1)]j ° lb * = ^ {21 + 1)[a 



(3) 
(4) 



The constant uj^ introduced above is the volume of the unit sphere in k-dimension. Wu = 



in 2D we will only need w = 1, uii = 2 and ui2 = 71". 



r(fe/2+i) 

The lowest order Hermite polynomials Hk(v) are listed below. The expression consists of two distinct contributions. The part which do not depend on 
the three different skewness parameters , S^ 1 ' , S^ 2 ' and signifies the MFs for a Gaussian random field. The other contribution 8Vk{v) represents the 
departure from the Gaunssian statistics and depends on the generalized skewness parameters defined in Eq.([6j and EqilOK 
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H-iW = y|expfyjerfcf-^J; H {v) = 1, H^u) = v, 
H 2 (u) = v 2 -l, H 3 (v) = u A -Zv, H 4 (u) = v 4 - Qu 2 + 3. 

ff.M -(-!>•«*(£) (5) 

Various moments Oj that appear in Eq.© can be expressed in terms of the power spectra Ci and the observational beam bi. The moment ao is a spacial 
case which correspond to the variance. The quantities o\, 02 are natural generalizations of variance and they put more weights on higher harmonics. The 
variance that will mostly be used are 0% = (\& 2 ) and a\ = ((V*) 2 }. 

The real space expressions for the triplets of skewness are given below. These are natural generalizations of the ordinary skewness S° that is 
used in many cosmological studies. They all are cubic statistics but are constructed from different cubic combinations. 



e(o) = _ <*!> „( 1)= J5^_J(*V$) (2) _ s(v^-v^v^) o ((W).(VV&)(V 2 M>)) 

*3 " 4 olo\ 4 a*o* a* °t ^ 

The expressions in the harmonic domain are more useful in the context of CMB studies where we will be recovering them from masked sky using 
analytical tools that are commonly used for power spectrum analysis. The skewness parameter S' 1 ' is constructed from the product field [\& 2 ] and [V 2 ^, 
where as skewness parameter S 1 ' 2 ' on the other hand relies on construction of [VP ■ V\&] and [V 2 \&]. By construction the skewness parameter S^ 2 ' has 
the highest weights on high / modes and <5r ' has the lowest weig hts on high I modes. The expressions in terms of the bispectrum B^^^ (see Eq. dl lb 
for definition) take the following form (see e.g. Hikag eetalJ I2OO8I) ): 



5 ( * ) — — B h i 2 i 3 Ii 1 i 2 i 3 bi 1 bi 2 bi 3 (7) 

h 

5' ( * 2v2 * ) = [ ll( ~ h + h{l2 + V + + l )] B hi2iJhhhKb h bi 3 (8) 

h 

^(Vf.wv 2 *) = _L [\h(h + 1) + l 2 (l 2 + 1) - i 8 (I 3 + l)]i 3 (i 3 + 1) + cyc.pcrm.] B hl2h I hl2h b h bi 2 b h (9) 



l hhh 



h 

h + l) (21-s + l) ( h h h \ , m 
4tt I J ' 1 ' 



We will use the same notation 6i 1 ! 2 ; 3 for the reduced bispectrum and the beam 6;. The number of subscripts however is enough to distinguish them. 

The bispectrum Bi x i 2 i 3 used here defines the three-point correlation function in the harmonic domain. In general a reduced bispectrum bi t i 2 i 3 is 
commonly used in certain situations which can directly be linked to the em flat-sky expressions. 



{ai imi ai 2m2 ai 3m , 3 ) c = [ m ^ m ^ j B h i 2 i 3 ; B h i 2 i 3 = Ii 1 i 2 i 3 bi 1 i 2 i 3 (11) 

The expressions for the MFs in Eq.© that we have discussed, depend on the one-p oint cumulants . Howev er it is possible to define power spectra 
associated with each of these skewness parameters following a procedure developed in Muns hi & Heavens! (2010). This will mean we can also associate 
a power spectrum with vj: which will generalize the concept of MFs in a scale dependent way. The power spectrum that we associate with MFs will 
have the same correspondence with various skew-spectra S\ that the MFs have with one-point cumulants or . The power spectra such defined will 
however have more power to distinguish various models of non-Gaussianity. This is one of the main motivation behind generalizing the concept of MFs, 
which is a number, to a power spectrum. 

The series expansion for the MFs can be extended beyond the level of the bispectrum. The next to the leading order corrections terms are related 
to trispectra of original fields and various derivatives constructed from the original field by using differential operations such as [V • V], [V 2 ]. These 
corrections are however expected to be sub-dominant in the context of CMB studies for the entire range of angular scales being probed. 

The results here correspond to temperature maps, which is a spin-0 object. It is possible to extend these results to the case of polarization analysis 
i.e. for spin-2 fields. Such results will also be context of weak lensing shear and flexions. A detailed analysis will be presented elsewhere. 

3 THE TRIPLETS OF SKEW-SPECTRA AND LOWEST ORDER CORRECTIONS TO GAUSSIAN MFS 

The skew-spectra are cubic statistics that are constructed by cross-correlating two different fields. One of the fields used is a composite field typically 
a product of two maps either in its original form or constructed by means of relevant differential operations. The second field will typically be a single 
field but may be constructed by applying various differential operators. All of these three skewness parameters contribute to the three MFs that we will 

consider in 2D. 

The first of the skew-spectra was studied by (Cooray 2001) and later bv lMunshi & Heavens! j2010 ") and is related to commonly used skewness. The 
skewness in this case is constructed by cross-correlating the squared map ^ 2 (Ci) with the original map ^ (Ci) . The second skew-spectrum is constructed by 
cross-correlating the squared map [^ 2 (f2)] against [V 2 >l/(f2)]. Analogously the third skew-spectrum represents the cross-spectra that can be constructed 
using • V#(f2)] and [V 2 *(n)] maps. 
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Figure 3. The skew-spectra sj~ (left-panel), (middle-panel) and S 1 , (right-panel) are plotted for various smoothing beam (see Eq.(l6j)as function of the harmonics I. 
The model for the primordial non-Gaussianity is assumed to be a local model for these plots Eq. j441 .The underlying cosmology is assumed to be that of WMAP7. The three 
FWHM for the three beams considered are 9 S = 10' (solid), 25'(dot-dashed), 50' (small-dashed), lOO'(long-dashed). The parameter fjST describing the normalization 



of local model of non-Gaussianity is set to unity / 
S* = £,(2i + l)S*. 
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Figure 4. Same as previous figure but for the equilateral model of primordial non-Gaussianity as defined in Eq. fl45t . The choice of smoothing angular scales (beam) are 
same as previous figure. The normalization parameter for equilateral bispectrum fWT % is set to unity /S^ 1 = 1 ■ We consider three different beams as indicated 8 S = 10' 
(solid), 25'(dot-dashed), 50' (small-dashed), lOO'(long-dashed). The shape of the two skew-spectra as derived from equilateral and local model are different and can be 
used to differentiate their contributions in real data. 
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These set of equations constitute one of the main results in this paper. The matrices here denote Wigner-3j (3j) symbols, bi represents the experimental 
beam bi = exp{— 1(1 + l)/2cr 2 } withcrj, = FWHM/V8 In 2 for a Gaussian beam. Each of these spectra probes the same bispectra Bu 1 i 2 with different 
weights for individual triplets of modes that specifies the bispectrum (7, hh)- These triplets of modes specifies a triangle in the harmonic domain. The 
skew-spectra sums over all possible configuration of the bispectrum keeping one of its side I fixed. For each individual choice of / we can compute the 
skew-spectra S, . The extraction of skew-spectra from data is relatively straight forward. It consists of construction the relevant maps in real space either 
by algebraic or differential operation and then cross-correlating them in the multipole space. The issues related to mask and noise will be dealt with in later 
sections. We will show that even in the presence of mask the computed skew spectra can be inverted to give a unbiased estimate of all-sky skew-spectra. 
The presence of noise will only affect the scatter. We have explicitly displayed the experimental beam bi in all our expressions. 

To derive the above expressions, we first express the spherical harmonic expansion of fields the [V 2 ^/(f2)], [V*l/(f2) ■ V*l/(f2)] and [^ 2 (f2)] in terms 
of the harmonics of the original fields ^i m . These expressions involve the 3j functions as well as factors that depend on various U dependent weight 
factors. 

[V 2 *(S7)] im = f dhYC m {Cl) [V 2 *(^)] =-J(I+ l)*, m 
[* 2 (^)] im = f d£lY t * m ((l) [* 2 (fi)] = ^(-l) m f hmi *i 2m J hl 

[V*(fi) • Vtf(A)]i m = /'dfty^(ft)[V*(ft)-V*(A)] = 53*i imi *i 2 m a y dClYCaitl) [VY hm2 (Q)-VY l2m2 {Cl)} (17) 

lirai 

= \ y £,[h(ix + i) + h{h + ±)-i{i + i)\ J d6Y^n)Y; imi (n)Yi 2m2 (fi) 

We can define the power spectrum associated with the MFs through the following third order expression: 



h h I 

mi m2 — m 



Vi z) = 5>*M 2 ' + 1) = ~ $>I + 1) {fflM + jS™H( V ) + S?>H k .,(,) + ■■■}. (19) 

The three skewness parameters defines the triplets of Minkowski Functional. At the level of two-point statistics, in harmonic domain we have three 

(3) 

power-spectra associated with Minkowski-Functional Vt that depend on the three skew-spectra we have defined. We will show later in this paper that 

the fourth order correction terms too have a similar form with an additional monopole contribution that can be computed from the lower order one-point 

terms such as the three skewness defined here. The result presented here is important and implies that we can study the contributions to each of the MFs 

Vk {y) as a function of harmonic mode I. This is a especially significant result as various forms of non-Gaussianity will have different I dependence and 

hence they can potentially be distinguished. The ordinary MFs adds contributions from individual I modes and hence have less power in differentiating 

various contributing sources of non-Gaussianity. This is one of main motivations to extend the concept of MFs (single numbers) to one-dimensional 

objects similar to the power spectrum. 

In Figure-[2]we have plotted the variance parameters [<r 2 ]; and [a 2 ]; for various smoothing beam (assumed Gaussian). The four different FWHM 

that are being considered are 9 S = 10', 25', 50' and 100' respectively. These parameters only depend on the underlying CMB power spectra and the 

beam. They are used as a normalization parameters while constructing the MFs from the generalized skewness parameters (see Eq.UOIl). The specific 

model that was considered for this particular plot is called local model (defined in Eq.J44t). The skew-spectra that correspond to the other popular model 

for primordial bispectrum is known as equilateral model (see Eq.l|45ll). We plot the corresponding skew-spectra in Figure-|4]as a function of angular 

harmonics I. The shape of the skew-spectra for various model can provide the necessary clues to the origin of primordial bispectra while analyzing CMB 
mapsin Figure-[3] we have presented the three different skew-spectra defined here Si as a function of the harmonics I. The skew-spectra for a generic 

bispectrum is defined in Eq.J 12t. Eq.lll3t and in Eq.lll5t. The skew spectra are sensitive to the smoothing bi moreover the skew-spectra at a given I depend 

on the bispectrum defined over the entire range of I being probed. The normalization parameters for these plots are set to be equal to unity; i.e. we take 

/nl ~ 1 an d /jvzT = !• The skew-spectra will scale linearly with these parameters. 

The point sources act as a contamination in estimation of the primordial bispectrum. The skew-spectra for the point source are plotted in Figure-[5] The 
specific model for the bispectrum that we consider is given in Eq.l l5U . The point sources are expected to dominate at higher I values. The normalization 
for the point source bispectrum is set by the parameter bps- In our study we take bps = 10 -27 /iK s 

Next, we consider the higher order corrections to the MFs. These corrections takes contribution from the trispectrum. Corrections to the individual 
MFs can be expressed in terms of a set of four kurtosis which are formed from the trispectra. These kurtosis are one-point estimators and they differ 
the way sample the individual modes of the trispectra defined by the quadruplet of harmonic number {h}; i = 1,2, 3, 4. These generalized kurtosis 
parameters (denoted by K l ) can be generalized to kurt-spectra, denoted as K\, in a manner very similar to the skew-spectra. These kurt-spectra can be 
used to express the next order corrections to the power spectra associated with MFs. 



4 THE QUADRUPLET OF KURT-SPECTRA AND NEXT TO LEADING ORDER CORRECTIONS 

In a perturbative analysis, the leading order terms that signify non-Gaussianity in the analysis of MFs depend on the bispectrum or equivalently a set 
of skewness. The next to the leading order order correction terms depend on a set of kurtosis parameters K^' that are fourth order statistics and are 
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Figure 5. Same as previous figures but for point sources. For the point sources we have taken bps = 10 _27 /iK" 3 . The panels represent the three different skew-spectra 
S^ ' (left-panel), (middle panel) and Sp' (right-panel) respectively. The skew-spectra are sensitive to the resolution of the maps being analyzed as they are integrated 
measures. This is related to the fact that it's value at a specific harmonics I takes contribution from all possible modes, i.e. the entire range of I values. The point source 
contributions dominate the primordial non-Gaussianity at high I. 



analogs of the skewness parameters S^' we ha ve defined above. In general the kurtosis parameters are collapsed fourth order one-point cumulants and 
probe the trispectrum with varying weights (see Muns hTet alj d201(j|) for a more detailed discussion on fourth order one-point cumulants, their two-point 
counterparts or the cumulant correlators and related harmonic space description). The four different kurtosis parameters that are related to the MFs are a 
natural generalisation of the ordinary kurtosis K^ ' which is routinely applied in many cosmological studies. We will denote these generalised kurtosis 
parameters by K^ 1 '; i = 1, 2, 3. These parameters are constructed from the derivative field of the original map map ^(Cl) and its derivatives [V*(S7)] 
and [V 2 *(0)]. 
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The subscript c correspond to the connected components which indicates that all Gaussian unconnected contributions are subtracted out, these include 
both noise as well as the signal contribution. The eavaluation of these moments are relatively easy in real space for a pixelised map. It involves taking 
derivatives of beam smoothed maps. The corresponding power spectra associated with these fourth order moments are constructed by cross-correlating 
appropiate maps in the harmonic domain and are easy to implement numerically (Eq.d40b provides exact expressions for the corresponding power spectra 
or the kurt-spectra.). 

The next to leading order corrections to the MFs involve these K^s as well as the product of two skewness parame ters (i.e. the terms such as 
[5 (0) ] 2 , [S (0) S (1) ] or [S (1) S (2) ]) defined previously in the context of leading order non-Gaussian terms ( lMatsubarj|2O10h . 
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The analytical modelling of four-point correlation functions is most naturally done in the harmonic domain. They are described by the angular trispectrum 
T lsu( L ) wmc h is defined through the relation (ai imi ai 2m2 ai 3 m 3 ai im4 ) e = Y,l ^a^iaU^^U The trispectrum T^ 1 / 2 (I) is expressed in terms 
of the reduced trispectrum P l \ 2 (I). Following expression was introduced bv lHul l l200lh and encodes all possible inherent symmetries. 



Tti:(i)=pi£(i)+{2i+i) 
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\h+h 



h h I 
h la I 
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h h I 

la h I 



(25) 



The matrices in curly brackets represent 6j symbols which are defined using 3j symbols (see lEdmondsl d 19681) for more detailed discussions). The entites 



P/ 1 / 2 (I) can be further decomposed in terms of the reduced trispectrum r/ 1 ;' 2 (I): 



Pllilil) = rll\l{l) + {,-lf u rll\l{l) + {-lf L r\lll{l) + (-1) 



M0; 



h+h + L; 



h + h + L. 



(26) 
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Each individual model for primordial non-Gaussianity makes a specific prediction for the reduced trispectrum which can be used as a fingerprint to rule 
out many possibilites. 

Next we will introduce three additional trispectra that are constructed using different weights to the original trispectra TJ^^L) and differ the way 
they weight various modes,that are specified by a particular choice of the quadruplet {k}. 

[T^tXiL) = [T w ]\^(L) = 1 [h(h + 1) + h(h + 1) + h(h + 1) + k(U + 1)] Z££(L); (27) 

l Ti2) ]\lu( L ) = i I L ( L + !) ~ ('id + !) + ^ + iJX'sCs + 1) + + !))] ^4 ( L ); C 28 ) 

P^]!^) = \ KHh + 1) + i 2 (i a + 1) - i(i + l))(is(i3 + 1) + Utf* + 1) - L(L + 1))] t££(L). (29) 

The four generalised kurtosis and the related kurt-spectra we have defined above can now be expressed in terms of these generalised tri-spectra as 
follows: 

tf^EE^aU^WW; =Y^^\\\l{l)Ji^iJi 3 i^ ^)(2I + l)Jf«=Jf«. (30) 

The fourth order expressions for the power spectra associated with MFs v\ can be obtained by replacing all one-point K^s with their two-point 
counterparts in Eq.J24t, The contributions from the skewness parameters only contribute in the monopole terms. The extraction of the kurt-spectra 
in the presence of a mask can be carried out analogus to the skew-spectra. 

On a different note, the power spectra associated with the kurtosis or kurt-spectra are discussed in detail for a scalar field dMunshi et alfeoiol) . two 
different types of kurt-spectra were introduced in the context of analysis of CMB Temperature maps. These two kurt-spectra K\ ' and K\ ' both 

(2 2) 

samples the relevant trispectrum. The first of these is constructed from cross-correlating the squared map with itself and is denoted as K\ ' . The other 

(3 1) 

kurt-spectra is constructed from cross-correlating a cubic map with the original map K \ .In general different set of maps can also be used to form 

(2 2) 

squared and cubic combinations which will probe a mixed trispectra. In the present context we are interested in the spectra K\ ' as the construction of 
K) ' will involve gradient maps and are more complicated to analyse in a coordinate independent way. 

The physical meaning of these kurt-spectra can be understood more easily in the harmonic domain. As mentioned, each individual modes of the 
trispectra is characterized by a spcific choice of set of modes U that defines it. These modes each constitute the sides of a quadrangle whose diagonal is 
specified by the quantum number I. The kurt-spectra that we have considered here take contributions from all possible configurations of the quadrangle 
while keeping its diagonal / fixed. 

The estimation of the kurt-spectra from a real data is relatively easy and follows the same methodology as that of the skew-spectra. The first of 
these kurt-spectra is extracted by cross-correlating the squared field ^ 2 (Cl) with itself. The spectra Jf' 1 ' is constructed by cross-correlating Vl/ 2 (f2) 
against $V 2 $. The other two kurt-spectra can likewise be constructed. In each such construction a scalar map from a product field is generated before 
it cross-correlated with another such map. The reason we have not constructed the other kurt-spectra namely K\ ' is simply realted to the fact that 
scalar maps are not possible for all such cubic combinations e.g. those involving gardient fields i.e. in a coordinate independent way. However such 
constructions are indeed possible using the spinorial formalism. 

The corrections to the power spectrum associated with the MFs now can be written in terms of the i.e. the kurt-spectra and the various Hermite 
polynomials as introduced above. The contributions that higher powers or lower order statistics such as skewness will only contribute to the monopole 
term for every MFs. However the higher multiploes will involve contribution from various Kurt-spectra as indicated. Different specific choice of Trispectra 
will therefore lead us to completely different power spectra associated with the MFs and can help to distinguish various models of non-Gaussianity. 



5 ESTIMATORS AND THEIR COVARIANCE 

The results derived above correspond to the all-sky and no-noise situation. However, in reality often we have to deal with issues that are related to 
the presen ce of a mask and ( inhomogeneous) noise. To correct for the effect of a ma sk and the noise w e will follows the Pseudo-CV (PCL) method 
devised bv lHivon et al.l J200lh for power spectrum analysis and was later developed by dMunshi et al.l l2009) for analyzing the skew spectra as well as the 
kurt-spectrum dMunshi et aDl2010h . 

The partial sky coverage introduces mode-mode coupling in the harmonic domain. Individual masked harmonics become linear combinations of 
all-sky harmonics. The coefficients for this linear transformation depends on specific choice of masks through its harmonic coefficients. We will devise a 
method that can be used to correct for the mode-mode coupling . Our method is based on PCL approach which is often used in reconstructing unbiased 
power spectrum. If we have a generic field A(Q) and B(Cl) we denote their harmonic decomposition in the presence of mask w(fl) as Ai m and Bi m . 
Notice that the mask is completely general and our results do not depend on any specific symmetry requirements such as the azimuthal symmetry. The 
fields A and B may correspond to any of the fields we have considered above. In a generic situation A and B will denote composite fields and the 
harmonics Ai m and Bi m will correspond to any of the harmonics listed in Eg ,d 1 8ft. 



.4 
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= j dClYi m (Cl) w{Cl) A{Cl); Bim = J d£lYC m {(l) w(£l)B{(l) 
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The above expression relates the masked harmonics denoted by Ai m and Bi m with their all-sky counterparts Ai m and Bi m respectively. In their derivation 
we use the Gaunt integral to express the overlap integrals involving three spherical harmonics in terms of the 3 j symbols (Edmonds 1968). The expressions 
also depend on the harmonics of the mask Wi m - If we now denote the (cross) power spectrum constructed from the masked harmonics and denote it by 
Si and its all-sky counterpart by Si we can write: 



qA,B _ 

bl ~ 21 



l^i^SL; Sf' B =Y / MwSf' B ; M„, = £ I^ vt \vn»f\ {A, B} e {*, * 2 , (V* ■ V*), V 2 *} . (33) 



In above derivation we have used the orthogonality properties of the 3j symbols. It is interesting to notice that the convolved power spectrum estimated 
from the masked sky is a linear combination of all-sky spectra and depends only on the power spectra of the mask used. The linear transform is encoded 
in the mode-mode coupling matrix M U i which is constructed from the knowledge of the power spectrum of the mask. In certain situations where the 
sky coverage is low the direct inversion of the mode mixing matrix M may not be possible sue to its singularity and binning may be essential, based 
on these results it is possible to define an unbiased estimator that we denote by S A ' B . The noise due to its Gaussian nature, do not contribute in these 
estimators which is remain unbiased. However presence of noise is felt due to increase in the scatter or covariance of these estimator which can be 
computed analytically: 

S?' B = £[M-V^' B ; {5S?< B 8S-^ B ) = Y, M^iSS-^SSt^M-),; {S A ' B ) = Sf' B ; {A, B} g {*, * 2 , (V* ■ V*), V 2 *} . 

I' LL' 

(34) 

Notice that the mode coupling matrix M is independent of the particular choice of the skew-spectrum that we are interested. Hence the same coupling 
matrix can be used to extract the power spectrum associated with the MFs as the MFs are constructed from the linear combinations of generalized 
skew- spectra. 

I' LL' 

In the perturbative limit where the MFs are computed using the analytical results in Eq.«(4j the information content of MFs is same as that of the 
bispectrum. The three skewness parameters that were defined in Eq.(|6) all probe the same bispectrum with varying weights. 
The variance {8S A ' B 8S A ' B ) of various estimators can be constructed using the following procedure: 

(SS A ' B 5S A ' B ) = C A ' A C B ' B + [S A ' B ] 2 (36) 

c v-v,v.v = ^CCAHh + 1) + Hh + 1) - 1(1 + I)] 2 /?,,-; Cr 2 '* 2] = Y, C l' C i»&l"-> Cr a * ,va * ] = l 2 (l + 1) 2 C (37) 
vi" vi" 

The symbols C A ' A denote the power spectrum of the fields A(Cl); A(&) is a generic field that are used for the construction of generalized skew-spectra. 
The derivation of the covariance depends on a Gaussian approximation i.e. we ignore higher order non-Gaussianity in the fields. Ci is the ordinary CMB 
power spectra it also includes the effect of instrumental noise Ci = Cf + Cf . For a survey with homogeneous noise, ignoring the effect of beam we can 
write C; = Q, p a% where Q p is the pixel area and ajv is the noise variance. In a regime when noise contributions dominate the MFs can be approximated 
by a Gaussian approximation. 

The estimators for various skew-spectra are expected to be correlated to a certain extent. These can be expressed using following expression: 

(SS Al ' Bl SS A2 ' B2 ) =C Al ' A2 C Bl ' B2 + S Al ' B2 S A2 ' Bl ; {A 1 ,B 1 ,A 2 ,B 2 } g {*, * 2 , (V* • V*), V 2 *} (38) 

The above results are sufficient to compute the lowest order corrections to MFs due to the presence of non-Gaussianity as well as the scatter in the 
estimates in the presence of realistic mask and noise. 

However the above method is also applicable to compute the leading order corrections to these terms, i.e. the kurt-spectra. These can be estimated 
using the results in Eq.ll22t. e.g. to construct the first of these kurt-spectra we have to cross-correlate the squared field ty 2 with it itself. Other kurt-spectra 
are similarly constructed by cross-correlating quadratic constructs which involve the derivative fields (see Eq.d29t and Eq.l|30|l): 

m m 

Kfa = 27TlE[* 2Uv2 *]*™; ^ 2i, = 2[rrE[*V[*V 2 *]? m ; Kf = - [ l I ^t| V *| 2 ] im [|V*| 2 ]L- (40) 

mm m 

We have defined two different estimators K 2a and Kf a that can be jointly used to construct the kurt-spectra Kf. To compute the Gaussian contributions 
equivalent estimators needs to be define that works with simulated Gaussian maps with same power spectrum. The resulting Gaussian estimates will have 
to be subtracted from the above estimates. Equivalently, the Gaussian compone nts of the kurt- spectra can be computed from Eq.l|29t by replacing the 
generic trispectra with its corresponding Gaussian counterpart G\ lt2 (L) (see e.g. lHul (|2000) and lHul l l200lh ): 

G\±(L) = (-1)'^ y/(2h + 1)(2Z 3 + rjC h C h SL S hh 6 h i 3 + (2L + l)C h C h [{~l) l2+l3+L Si^i, + 6 hu 5, 2 i 3 ] . (41) 
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where Ci is the power spectrum including noise. 

The treatment for the masked sky will follow in an exactly same manner. For individual the unbiased estimators can be recovered using exactly 
the same mode-coupling matrix M«/ introduced before in Eq.<l34b: we have for the masked kurt-spectra — M^, 1 . The auto- and the covariance 
of these estimators can also be estimated using obvious generalization of Eq.f34t and Eq.l|38t respectively. However previous studies have concluded that 
for realistic values of /jvl the corrections terms that involve the kurt-spectra will be almost negligible. 



6 MINKOWSKI FUNCTIONALS AND THE CMB SKY 

The discussion so far has been completely generic and is applicable to arbitrary random 2D field on the surface of the sky. We will specialize the 
discussion in this section to the case of CMB sky. CMB observations are the cleanest probes of primordial non-Gaussianity. The angular multispectra for 
the temperature fluctuations sample the 3D multispectra of the inflationary potential. 

Given a specific form for the primordial non-Gaussianity It is possible to compute the MFs for the observed temperature perturbations. The non- 
Gaussianity in the CMB sky can be a direct manifestation of the non-Gaussianity in the seed perturbations generated during the inflation. The non- 
Gaussianity in inflationary potential is most easily characterized in the Fourier domain. The decomposition provides 3D Fourier coefficients $(k) for 
the wave vector k. The following expression links the curvature fluctuations in Fourier space generated during inflation < l > (k) with spherical harmonic 
coeffi cients of temperature anisotropy a; m with the help of radiation transfer function A;(fc) for the temperature fluctuations dWang & Kamionkowskil 
2000). The angular power spectrum for the temperature fluctuations Ci = (a; m a* m ) can be expressed in terms of the power spectrum of the 3D 
perturbations in the potential field (<E>(ki)<£>(k2)} = (27r) 3 5 3 D(ki + k 2 )P$(fci). 



d 3 k 



a lm =M-iYJ j-^9(k)At(k)Y tm {k); C t = - J k^dkP^{k)Af{k). (42) 

A Gaussian sky can be described statistically just by its angular power spectrum Ci. The lowest order departure from the Gaussianity is described 
in terms of the bispectrum angular bispectrum. The general form for the 3D bispectrum for the inflationary potential $ is given as B<s>{ki, fei, k 3 ) = 
(27r) 3 (53_D(ki +k2 + ks) F${ki, &2, ^3). In general translational invariance enforces the momentum conservation in the Fourier domain will introduce the 
factor 5{id (ki + k 2 + ks) involving the 3D Dirac delta function 830- The kernel F®{k\, fa,k 3 ) therefore is the amplitude of the bispectrum associated 
with each triangular configurations involving the wave vectors ki. Various early Universe scenarios differ, the way they specify F&(fa, k 2 ,k 3 ). The 
reduced angular bispectrum defined in Eg. ill It can be expressed in terms of F$(ki, fa, £3): 



^\ J drr 2 j k 2 dkij h (fcir)A ;i (fcir) j k 2 dk 2 j h {far) A h {k 2 r) j k%dfaji 3 {far) A h (fc 3 r)P $ (fci, ki, fa). 



(43) 



Models of inflation can largely be divided into two different categories. The first class of models is known as the local mod eld45t 1 K omatsu & SpergeJ 
l200ll ; lMedeiros & Contaldj200d : ICreminellil2003l : ICreminelli et alj 20061 : ICabella et al. l2006r . lLiguori et al.l200llSmifh. Senatore & Zaldarriagj2009l) . 
In these models the contrib ution to the bispectrum is maximum for the squeezed confi gurations i.e. when k\ ;§> fa, fa. The other main class of models are 
called equilateral models dChen. Huang & Kachru 2006; Ch en. Easther & Limll2007r) . In these class of models the maximum contribution correspond to 
a configuration where all wave vectors have similar magnitudes ki ~ k 2 ~ fa i.e. a equilateral configuration of the triangle representing the bispectrum 
in the Fourier space. It is important to note that unlike the loca l model the equilatera l model can not represented by product of separable functions. 
However approximate separable forms do exist in the literature dCreminelli et alj|2006l: ISmith & Zaldarriagall2006l) . We quote the results for the CMB 



non-Gaussianity that arise in the context of local or equilateral models iKomatsu & Spergel BoOlh 



bi°J 2h = 2/ivL / rdr[/3 h (r)fr 2 {r)a h (r) + cyc.perm] (44) 



KZh = 6 fNL J rdr[-a h {r)(3 h {r)J3 h (r) - 26 h {r)6 h {r)S h (r) + p h {r) 7h (r)S h (r) + cyc.perm.] . (45) 
The following functions are useful in analytical expression for the bispectrum and trispectrum dCreminelli et al.l l2006): 



ai{r) = - [ k 2 dkP^{k)j l {kr)A l {k); ft(r) = - f fc 2 dfcP 4 >(fc)j ; (fcr)A ! (fc); (46) 

7T J TV J 

Mr) = k 2 dkPl /Z {k)j l {kr)A l {k); 5i{r) = k 2 dkP 2/3 {k)j t {kr)A t {k); (47) 

F L {n,r 2 ) = - [ k 2 dkPl /Z {k)j l {kr 1 )ji{kr 2 ). (48) 



Here ji is spherical Bessel function, $ is the inflationary potential A; (fc) is the radiation transfer function which can be computed using the publicly 
available software such as cameB or CMBFASllE In addition to the bispectra the re duced CMB trispectrum in the local model can be expressed in 
terms of these functions as (Hu 2000. l200ll ; lKomatsu & SpergefeOOlklKogo et aljioO^) : 



3 http://camb.info/ 

4 http://www.cmbfast.org/ 
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T illl( L )= 4f§j h hi x i 2L hi 3 i iL J rfdn J r2dr 2 F L (r 1 ,r2)a h (r 1 )f3 l2 (r 1 )a h (r2)l3 u {r 2 ); 

+gNhhi lhL hi 3 i A L J r 2 dr/3 i2 (r)/?i 4 (r)[^i 1 (r)/3 i3 (r) + ^ i3 (r)ft 1 (r)]. (49) 

For detailed discussions about various issue s related to the symmetries and modeling of the CMB trispectrum see fau & Okamo"toll2002l ; iHufcOOd 1200 ll ; 
iKomatsu & SpergefcOOlklKogo et alj2006l) . 

In the limit of low multipole the perturbations are dominated by Sachs-Wolfe effect. In this regime the transfer function A; (k) takes a simpler form; 
A(fc) = 1/3 ji (fcr*). Here r» denotes the time elapsed since cosmic recombination = t)q — r jdec- The specific form for the reduced bispectrum bj-, and 
reduced Tj*/' (L) are completely specified by the angular power spectra C; in this limit (Hui 200(il. l200ll ; lKomatsu &~S pergel 200ll; lKogo et alj2006h : 

bfZh = -6fNL(C h C h +C h C l3 +C l3 C h ); [r 1 ^ (L)] sw = 90^(4 f 2 m C L + g NL (C h + C h ))I hl2L I hUL . (50) 

The MFs in this limiting case will therefore be completely specified by the angular power spectra. The fact that MFs can be specified completely by the 
bispectrum for arbitrary smoothing scales simply means that for local as well as equilateral models the parameter /jvl is sufficient to fully characterize 
the MFs. 

The point sources too contribute to the CMB bispectrum. The point source bispectrum is given by: 



^1*2*3 = bps. (51) 

The exact value of the amplitude bps depends on the limiting flux used in a specific survey. In our study we have taken bps = 10 -27 fiK 3 . The 
results of our computations are plotted in Figure-l[5]l. We have considered three different Gaussian beams as indicated 9 S = 10', 25', 50', and 100'. At 
higher resolution I the point-source contribution will dominate primordial bispectrum. The skew-spectra Sf which puts more weights on smaller angular 
scales are more dominated by point source contributions. 



7 CONCLUSIONS 

Different approaches are essential while testing non-Gaussianity as they all exploit different statistical characteristics. There is no unique approach that 
can be adopted to describe or parametrize non-Gaussianity in a complete manner. Testing of non-Gaussianity therefore must be done using a battery of 
complementary techniques. Each of these techniques has a unique response to the real world issues such as the sky-coverage and instrumental noise. 
Any robust detection therefore will have to involve a simultaneous cross-validation of results obtain from independent methods. The most common 
characterization of non-Gaussianity involves a study of the bi-spectrum which represents lowest order departure from Gaussianity. Higher order non- 
Gaussianity can be studied using its higher order analogs i.e. the multi-spectrum. In contrast the topological estimators that we have studied here carry 
information of all-orders though in a collapsed ( one-point) form. Analytical results for MFs for a Gaussian field is well understood; which forms the 
basis of all non-Gaussianity studies iTomitj 1986). There has been several p revious studie s on ex traction of the MFs from the CMB data that rely either 
on simplification of radiative transfer using the Sachs-W olfe limit (see e.g . Hikage et al.l (|2006)) and using a perturbative approach based on a series 



expansion of the MFs that can be studied order by order dMatsubarall2003 , ll994l) . The MFs have also been studied using elaborate computer intensive 
non-Gaussian simulations (Kom atsu et al. I l2003l ; ISpergel et alj|2007h . Most of these studies were done using a specific model for the non-Gaussianity 
namely the local model introduced here in Eq.lJTJ which is parametrized by the parameter }nl- 

One of the main difficulties faced by one-point estimators, which also include the MF based approaches, is their inability to differentiate various 
sources of non-Gaussianity. The one-point estimators all effectively compress all available information to a single number. This has the advantage of 
increasing the signal to noise ratio but it loses the ability to differentiate various sources of non-Ga ussianity. In any study of pri mo rdial non-Gaus s ianity it 
is of the utmost importance to avoid any cross-contamination from secondary sources (see e.g. iGoldberg & Spergell dl999allbl) ; ICoorav & Hul feOOOh ). 
In recent studies invol ving MFs, certain levels of disagreement has been noticed with studies that use bi-spectrum to probe non-Gaussianity (see 
e.g lHikage et ail (2006)). Given the MFs based approaches only directly probe the bispectra, as the contributions from the higher order multi-spectra 
are sub-domina nt, it is important to unde rstand various reasons of these disagreements. 

Recently iMunshi & Heavensl d2010h have developed a new technique to study non-Gaussianity. Instead of one-point estimators, e.g. the skewness, 
their method relies on a spectrum or skew-sepctrum which is the Fourier transform of two-point objects known as cumulant correlators. These skew- 
spectra do not compress the all available information from the study of a bi-spectrum to a single number and their shape can help to distinguish among 
various sources of non-Gaussianity. Exploiting the perturbative expansion of the MFs it is possible to show that in the leading order of non-Gaussianity, the 
MFs depend on three generalized skewness parameters. We extend the concept of the skew-spectra to the study of these MFs and introduce a generalized 
skew-spectrum associated with each of these skewness. This allows us to introduce a series of numbers associated with each of the MFs. We call these 
objects the power spectrum associated with MFs. The advantage of cross-checking the contributions to MFs using the skew-spectra is they provide a 
method to test any contamination from sources which are simply cross-contaminants from secondaries or foregrounds. The two methods that compute 
the MFs should provide identical results. The method based on the skew-spectra are simpler to implement once the derivative fields are constructed. 
These methods are similar to moment based approaches for studying non-Gaussianity and hence can provide valuable basis for cross-comparison. We 
have shown that this can implemented in a model independent way. Our method is based on a Pseudo-C; approach and can handle arbitrary sky coverage 
and inhomogeneous noise distributions. The PCL approach is well understood in the context of power spectrum studies and its variance or scatter can 
be computed analytically. We provide generic analytical results for computation of scatter around the individual estimates. It is possible indeed to go 
beyond the lowest level in non-Gaussianity. However, it is expected that such correction will be sub dominant at least in the context of CMB data analysis. 
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Nevertheless we include the spectrum associated with the next order correction t erms that were intro duced by (Mats ubaralkoiOl) . These terms represent 
kurtosis and the corresponding spectra are known as kurt-spectra. In their study (Munshi et al. 2009) introduced two set of kurt-spectra. Adopting their 
method we show that generic two-to-two kurt-spectra can be extracted from the data without introducing any additional complication. To summarize: 

• In the context of CMB, the study of MFs is equivalent to the study of bispectrum. 

• Higher order terms which depend on trispectrum make negligible contributions. 

• However the MFs do not probe the full bispectrum but only weighted sums of modes or the three generalised skewness parameters. 

• We have defined three generalized skew-spectra associated with each of these skewness parameters. 

• Study of these skew-spectra can replace the study of MFs. 

• The skew-spectra can all be probed for arbitrary mask and noise. 

• The unbiased estimators can be constructed which can work in the presence of partial sky coverage. 

• Their variance can also be computed analytically there by avoiding the use of non-Gaussian simulations completely. 

• The MFs can be constructed from the knowledge of generalized skew-spectra and can be compared with the results from real space analysis. 

• Additional skew-spectra can be constructed using various derivative fields constructed from the original filed and can provide additional cross-checks 
or even help to separate various components of non-Gaussianity. 

• Generalized skew-spectra can be used to separate individual components of NGs. 

The approach we have developed here to study MFs in this paper will be most useful in case of weak lensing studies where real space analysis is difficult 
because of complex survey geometry. A detailed analysis will be presented elsewhere. 

The generic results derived here will also be applicable to other areas which include analysis of redshift surveys (3D) or projected surveys (2D) that 
include weak lensing surveys or SZ surveys. Results of detailed analysis will be presented elsewhere (Munshi et al. 2010 in preparation). It is however be 
important to keep in mind that unlike the other optimized estimators of Jnl based on bispectrum, the MFs are sub-optimal in nature. 

In this paper we have concentrated mainly on all-sky expressions that involve expansion of relevant fields using spherical harmonics. In cer- 
tain situations the flat-sky results may be more directly relevant. The flat sky approximations can be recovered by replacing the spherical harmonics 

Yim{9, <j)) ~ (— l) m y^ ^f-Jm[(l + exp(im4>). These results will be valid for high I ^> 1. Typically these will result in expressing the flat sky 

multi-spectra in terms of their all-sky counterparts i.e. C; ~ C(l + ^) for power spectrum and S; 1 ; 2 i 3 « Ii 1 i 2 i 3 B{h + \, h + \, h + §) for the bispectrum 
and an analogus relation holds at the level of trispectrum. 
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